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ABSTRACT 

We formulate and solve by semi-analytic means the axisymmetric equilibria of relativistic 
self-similar disks of infinitesimal vertical thickness. These disks are supported in the horizontal 
directions against their self-gravity by a combination of isothermal (two-dimensional) pressure and 
a flat rotation curve. The dragging of inertial frames restricts possible solutions to rotation speeds 
that are always less than 0.438 times the speed of light, a result first obtained by Lynden-Bell and 
Pineault in 1978 for a cold disk. We show that prograde circular orbits of massive test particles 
exist and are stable for all of our model disks, but retrograde circular orbits cannot be maintained 
with particle velocities less than the speed of light once the disk develops an ergoregion. We also 
compute photon trajectories, planar and non-planar, in the resulting spacetime, for disks with 
and without ergoregions. We find that all photon orbits, except for a set of measure zero, tend 
to be focused by the gravity of the flattened mass-energy distribution toward the plane of the 
disk. This result suggests that strongly relativistic, rapidly rotating, compact objects may have 
difficulty ejecting coUimated beams of matter or light along the rotation axes until the flows get 
well beyond the flattened parts of the relativistic mass distribution (which cannot happen in the 
self-similar models considered in this paper). 

Subject headings: general relativity: disks; general relativity: ergoregion; general relativity: frame 
dragging; general relativity: light bending 

1. Introduction 

It has been said that there are basically only two kinds of self-gravitating objects in astronomy, spheres 
and disks. It has also been said that general relativity is so beautiful, it has to be right. Thus, spheres 
and disks should be of as much natural interest to relativists as to Newtonian dynamicists. Yet because the 
beauty of general relativity comes at the steep price of great mathematical difflculty, for many decades the 
only known solutions to Einstein's field equations were ones possessing spherical symmetry. Any analytical 
attempts to study realistic rotating bodies relied on perturbation theory and the assumption of low angular 
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momentum. It took astrophysicists nearly fifty years after Einstein first wrote down the final form of his 
theory to find an asymptotically flat vacuum solution which has nontrivial angular momentum (Kerr 1963), 
and another thirty- three years to construct appropriate interior solutions (see, e.g., Lynden-Bell & Pichon 
1996). Some progress has been made in the study of relativistically rotating, axisymmetric objects with finite 
physical extension, but mostly in the limit when the disk is cold (Bardeen & Wagoner 1971, Lynden-Bell & 
Pineault 1978b), or when the material in it is taken in the form of two equal, coUisionless, counter-rotating 
sheets (Lynden-Bell & Pineault 1978a, Lemos 1989). The first case results in the mathematical simplification 
that the number of unknown metric functions reduces to three; the second, in the elimination of the dragging 
of inertial frames. 

However, the Newtonian analogs to cold disks are fraught with fierce dynamical instabilities (sec, e.g., 
Binney & Tremainc 1987). Such models cannot represent good approximations for realistic astrophysical 
systems (e.g., spiral galaxies, see Bertin & Lin 1996; or protoplanetary disks, see Adams & Lin 1993). 
Counter-rotating disks can avoid instability if they are sufficiently hot, but since such configurations have to 
arise as stellar-dynamical (coUisionless) rather than gas-dynamical (collisional) systems, relativistic analogs 
may have difficulty reaching the requisite degree of physical compactness. 

In the Newtonian studies of self-gravitatng disks, devices that have proven to have great mathematical 
utility are the assumptions of complete flattening and self-similarity (e.g., Mestel 1963, Zang 1976, Toomre 
1977, Shu et al. 2000). Razor-thin disks whose surface densities are power laws in radius r but which need 
not possess axial symmetry have solutions that can be found by analytical or semi-analytical means (i.e., 
involving nothing worse than the numerical integration of ordinary differential equations [ODEs]; see, e.g., 
Syer & Tremaine 1996, Galli et al. 2001). The gravitational collapse of such Newtonian models have elegant 
self-similar properties in spacetime (see, e.g., Li & Shu's [1996] study of the collapse of the axisymmetric 
singular isothermal disk). The relativistic analogs of such gravitational collapses, axisymmetric and non- 
axisymmetric, could lend valuable insight into issues of great contemporary interest in general relativity, 
such as the efficiency of gravitational radiation, or the possible formation of naked singularities. 

As preparatory work toward such applications, we wish to extend Lynden-Bell and Pineault's (1978b) 
work on cold, axisymmetric, relativistically rotating disks to obtain semi-analytical solutions of a family of 
relativistic disks parameterized by constant isothermal sound and rotation speeds, a and V. Self-similarity is 
then dimensionally still possible in the relativistic regime, because relativity introduces only another constant, 
the speed of light c, with the same units of velocity. Since no power of G, the universal gravitational constant, 
combined with c (or a or V) can yield a quantity with the dimensions of length (or time), it becomes natural 
and feasible to look for unbounded disk solutions where the surface density varies as a power of some 
appropriately chosen radial coordinate r. The mathematical consequences of this basic idea are developed in 
the paper as follows. Sections 2 and 3 derive the equations of axisymmetric stationary spacetime generated 
by a rotating disk. Section 4 describes our numerical strategy for the solution of the resulting ODEs. Section 
5 gives the results of the numerical integrations, and it also explores some properties of massive test particles 
placed in circular orbits in the disk plane. Section 6 considers the orbits of massless test particles (photons 
or neutrinos) in the general spacetimes of our models. Finally, in Section 7 we offer our conclusions and 
speculations. 
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2. Dimensionless Basic Equations 

2.1. Elementary Dimensional Considerations 

We characterize isothermal disks with flat rotation curves by two dimensionless parameters: the linear 
rotation velocity as a fraction of the speed of light, v = V/c, and the square of the isothermal sound speed 
as a fraction of c^, 7 = (a/c)'^. We nondimensionalize by adopting the unit of mass per unit length as c^/G 
= 1.35 X 10^^ g cm"-'^ = 0.677 Mq km~^. Notice that with c and G alone, we cannot define a characteristic 
mass per unit area (surface density), nor can we define a characteristic length. Thus, if r is a coordinate 
radius, with the physical units of length, we are naturally interested in disks with surface densities that are 
proportional to c^/Gr, i.e., with surface densities that are inversely proportional to one power of r. 

In general relativity it is possible to accomplish the equivalent nondimensionalization by working in the 

geometrical system of units where c = G = I. It is also possible, of course, to choose a radial coordinate tlp 
that does not have the units of length (see, e.g., Lynden-Bell & Pineault 1978b and Lemos 1989). However, 
to maintain self-similarity in the problem, tlp can be at best only some power (1/fc below) of r. Using tlp 
allows us to specify in advance that in going from pole to pole along a locus tlf = constant (for a slice 
at constant (f> and t), the associated polar angle 9 ranges from to tt, with the disk midplane located by 
symmetry at 6 = 7r/2. However, this adherence to normal convention comes only at the expense of making 
k a nonlinear eigenvalue of the problem. Using r and absorbing k into the definition of 6 eliminate the need 
for a complicated numerical procedure to find the value of k, but it puts the location of the midplane at a 
polar angle 6 = 60 ^ Tr/2. The former represents a considerable computational advantage, whereas the latter 
serves as a useful reminder that if we choose a coordinate r such that TZo^r represents the physical radial 
distance from the origin to the point (r, 9, (j)) in the disk midplane (with TZ having operational meaning as 
a proper radius because the origin is only mildly singular), then the distance along constant TZ (again for a 
slice at constant and t) from the midplane to the pole may not equal 7r7?./2 because of the distortion of 
the spatial geometry produced by the flattened mass distribution. 



2.2. The Metric 

Without losing generality, the metric that is stationary, axisymmetric and invariant under <j) ^ —cj) and 
t ^ —t may be written in geometrical units as 

= -e'^'^de + B2e-2-r.2 sin^ e{d(j) - udtf + e2(''-")(dr2 + r^dO^), (2-1) 

where ly, B, u) and \i are in general functions of r and 9. Since there is no fundamental unit of length, one 
might naively conclude that f, B, wr and ^ are functions of 9 only in a self-similar disk. This conclusion is 

premature and false. 

In the weak held limit, when the surface density and the rotation velocity are small, a cold Mestel disk 
of inflnite extent has the associated gravitational potential, 

(2-2) 

where D is a. fiducial length scale that contributes only an added constant to the potential and thus enters 
nowhere else in the problem. We will discard D in what follows. In the Newtonian limit, therefore, e'^ oc r", 
where n « when <^ 1 (with n having a different dependence on v and possibly also on 7 when 
the rotation and isothermal sound speeds in the disk are not small compared to c). This behavior - the 
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non-predetermined power-law of the gravitational distortion of time and the spatial geometry by the disk's 
self-gravity - is the source of the scaling relationships described by Lynden-Bell & Pineault (1978b) and 
Lemos (1989). According to the nomenclature of Barenblatt (1976), the situation is an example of self- 
similarity of the second kind. 

Equation (2-1) allows a spacetime that is self-similar under the transformation r ar, t a^~"'t. 
Then ds'^ — > a^ds^, and we can write the metric for < ^ < 20q as 

ds^ = -r^^e'^dt^ + r^e^^'-^'idcl, - r^-'e''-''Qdtf + e^-"" {dr^ + r^dO^), (2-3) 

where A^, P, Q and Z are functions of only. In the above metric, n is a pure number between and 1 that 
measures the depth of the disk's gravitational well. In particular, photons that emerge from the origin and 
reach some finite r have frequencies that are infinitely rcdshiftcd from their starting values. Self-similarity 
then implies that the same infinite redshift applies to photons that originate at any finite r and try to 
propagate to infinity (see §5). In a certain sense, therefore, the system constitutes an incipient black hole, 
one that will presumably acquire a growing point mass at the origin, with an accompanying horizon, if the 
disk is unstable and undergoes inside-out gravitational collapse (see Li & Shu 1996 for a Newtonian analog). 

It is instructive to compare our metric to the one used by Lynden-Bell & Pineault (1978b) and Lemos 
(1989): 

ds^ = -rl^e^'dt^ + r2|e2^-^(# - r™p-'=e^-^Qdt)2 + r^^^-^^ e^^^-"" {drip + rl^dx^). (2-4) 

The relationship between the two coordinate conventions can be made explicit by the following transforma- 
tion: 

r = r£p, e = kx, n=j, Z = Z^p~lnk. (2-5) 

While the two metrics arc completely equivalent, our choice turns out to be more advantageous in numerical 
implementation. For Lemos and Lynden-Bell & Pineault, the disk is located at x = '''/2, and m and k 
are eigenvalues of the problem. To find their values and the scaling of rotation velocity with density and 
pressure requires a three-dimensional shooting method, a nontrivial numerical task. For our metric, we have 
only one eigenvalue n. The other degree of freedom is embedded in the location of the disk, which we can 
find by satisfying certain jump conditions when we cross from the top to bottom surface. Later on, we will 
reparameterize the solution space to avoid even having to find n as an eigenvalue. 

Define the orthonormal tetrads for the locally non-rotating observer (Lemos 1989): 

^(0)" = (r-- r-'Qe^/^-P, 0, 0), (2-6) 

e(,)'' = (0,r-ie^/2-^,0,0), (2-7) 

e(,)'' = (0,0,e(^-^)/^0), (2-8) 

e(3;^(0,0,0,r-V^-^)/2). (2-9) 

The indices in parentheses label the basis vectors in {t,4>,r,9), and are raised and lowered with the fiat 
Minkowski metric r/ = diag(— 1, 1, 1, 1). The nontrivial Ricci components are 

2i?(o)(o)r'e^"'^ = ^99 + NePe + 2n(l + n) - {[(InQ)^ -P^ + Nef + {1- nf^ , (2-10) 



2i?(o)(i)r2e^-^ = Qee + Q9P9 - Q [Pee - Nee + {Pe - Nef + PeiPe - Ne) + 2(1 - n)] , (2-11) 
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(ii(o)(o) - i?(i)(i))r2e^-^ = Pee + + {n + l)^ (2-12) 
2ii(2)(2)r2e^-^ = Nee - Zee + Pe{Ne - Ze) + 2n(l - n) + ^^(1 - nf, (2-13) 
2i?(2)(3)r2e^-^ = (n + 1)^^ - 2nNe + Q^(l - n) [Pe - (InQ)^ - Ne] , (2-14) 

2 [i?(3)(3) + i?(o)(o) - i?(i)(i) - i?(2)(2)] r^e^-^ = 2Pe^e-iV2 + 4n2 + Q2 |[^^ + (i^Q)^ _ p,]2 _ _ i)2| . 

(2-15) 



2.3. Matter 

For the stationary thin disk, there is no radial or vertical motion. Hence we may write the four velocity 

(2-16) 

is the coordinate angular velocity and 



as 



where 



r-"e-^/2(l-t;2)-V2(i, 0,0,0), 



dt M* 



(2-17) 
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is the linear velocity of the fluid in the locally nonrotating frame. The physical significance of this quantity 
become clear when we project the four-velocity onto the locally nonrotating frame defined by the tetrad: 



= ( 



1 



Vi - ' Vi - v"^ 



,0,0), 



(2-19) 



i.e., the physics is that of special relativity in this frame, with v equal to a constant. In order to keep the 
disk infinitesimally thin, we will adopt the matter content of a fluid with vanishing vertical pressure. In the 
rest frame of the fluid, the stress energy tensor is given by 



t(a)(b) 



/£ 0\ 

P0 

Pr 

Vo 0/ 



d{9 - 0o). 



(2-20) 



Boosting into the tetrad frame, the nonvanishing components are (suppressing the argument of the delta 
function) 



T, 



mo) 



H0)(1) 



■(1)(1) 



1 2 ' 

{e+P4,)v 
p^ + ev"^ 



T(2)(2) =PrS- 

The Einstein field equations take the usual form 



-(a)(6) 



Via)ib)T 



(2-21) 

(2-22) 

(2-23) 
(2-24) 

(2-25) 
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where T = T^^'^^^ = {—e +P4, +Pr)S is the trace of the stress-energy tensor. 

The equation of motion on the disk is given by the vanishing divergence of the stress-energy tensor, 
T'^^.j^ = 0. For this particular configuration, only the /i = r component is nontrivial. The j^i = 6 component 
is proportional to the derivatives of the metric coefficients, which are discontinuous at the equatorial plane. 
By symmetry, this equation is identically satisfied once integrated through the plane. After some algebra, 
the r component reads 

n— — ^ 7 2 Qov{l-n)- 1 + (n + 2)pr + r— = (2-26 

1 — v' 1 — 1 — dr 

where Qo is the value of Q on the disk. With Qo, n, and v equal to constants, this equation implies that 
e, P0, and Pr all have the same radial dependence, which on dimensional grounds must be r~^, as can also 
be seen from their coupling to geometry through the Einstein field equations. Therefore the above equation 
becomes an algebraic one: 

Qov{l -n){e+ p^) + + ev^ - n{e + p^v^) - npr{l - v^) = 0. (2-27) 



3. Final Form of the Equations 

Define the reduced energy and stresses: 

£ = 87r^r2e^''-^°, p. = STr^r^e^""^", = g^^^r^e^""^", (3-1) 

1+n 1+n 1+n 

where a subscript denotes the value on the disk at 9 = Oq. Define a further rescaling, 

e = (1 + n)e, (3-2) 

and let a prime denote differentiation with respect to 0. Coupled to matter, the equation for P takes the 
form 

P" + P'^ + I = p,5{Q - Qo) (3-3) 

Away from the disk, we can solve this equation subject to the boLmdary condition e^(°-* = (so that a circle 
around the axis will have vanishingly small circumference as 0^0): 

P = In [sin 9] +C, P' = cot 6. (3-4) 

The constant C remains arbitrary, which enables us to set the boundary condition for other metric functions 
later. The solution (3-4) is only valid in the range < 9 < 9o where P is differentiable. For 9 > 9o, we 
can obtain the solution simply through symmetry considerations (i.e., the metric functions are even about 
the disk, while the derivatives are odd). Integrating equation (3-3) across 9 = 9o and combining the result 
with the second relation of equation (3-4), we get 

-2cot9o=Pr- (3-5) 



Let us confine our attention to < 9 < 9o. The rest of the field equations become 



N" + N'P' -h 



2n 
l + n 



{[{\nQ)' - P' + N'\ 



1 — n 
1 + n 



Pr + (e+P,^)Y3^ 



A, 



(3-6) 
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Q" + Q'P' - Q 



y ""l -1-N" + {N' - P'f - cot @N' 

[l + nY 



Qpr - 2{S + P4,)- ^ 

1 — 



A, 



Z' - Jl^N' + Q^L^ [P' _ (inQ)' - N'] = 0, 
1 + n 



Q^|[iV' + (lnQ)'-Pf -(i^ 



2P'Z' - iV' 



(l + n)2 



0. 



(3-7) 
(3-8) 

(3-9) 



where A = i5(6 — 9o) = S{6 — 6'o)/(l + n). Here we have dropped the equation that involves Z". It can be 
recovered by the equation of motion (2-27) owing to the contracted Bianchi identity. These equations can 
be simplified to give the dynamic equations, 



N" + N'P' + - {Q' - QP' + QN'f - Q2 

1 + n 



1-n 
1 + n 



(3-10) 



Q" + Q'P' - Q 



1 — n 
1 + n 



+ {N' - P'f - {Q' - QP' + QN' 



/x2 



= -{s:+P<l>) ^ A. 



and the constraint equation, 



1 - n 



4n 



, 1 — n 



4n^ 



(3-11) 



[QP' -Q'- QN']' - Q2 + T— - 20^' 7^ [Q^' -Q'- Q^'] - N + TT—^ = ^■ 

\l + nj 1 + n 1 + n (1 + n)^ 

(3-12) 

Now, Z has been completely decoupled from Q and N. We may use equation (3-12) to decouple Q from N 
as well, but this offers little advantage for numerical purposes. 



4. Numerical Implementation 

4.1. Boundary conditions 

We have already discussed the boundary condition for P. Unless r = 0, we expect space to be regular 
on the axis, which means the Ricci tensor remains finite there (and thus Riemann normal coordinates exist 
there). On the pole, P' = cotO diverges as 0~^. In order to have a regular solution for equation (3-12) 
there, we require 

N' = 0, = at e = 0. (4-1) 

The first requirement prevents the O geometry from having a cusp at the pole; the second discounts frame 
dragging on the rotation axis. One might naively expect Q' = on the pole as well from the term Q'P' 
in equation (3-11); however, this singularity is cancelled by QP' . In fact, we can use Q' on the axis to 
parameterize the solution space . With proper rescaling of t, r and e^, we can set A'' = on the pole (which 
means r" dt is the interval of proper time of an observer at 8 = 0). 

On the other hand, the delta functions on the right hand sides of the governing ODEs (3-10) and (3-11) 
signal a discontinuity in the derivatives of metric coefficients when we cross the plane of the disk. Similar 
to the method we used to find P, we integrate across the disk, and assume that the geometry is symmetric 
about the disk plane. Equations (3-10) and (3-11) then yield the following boundary conditions at 9 = ©q: 



Pr + {e+P,l>)- 2 



n' 1^- , ~ , Qo + Qov' + 2v 
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4.2. Scalar Two-Dimensional Pressure and Method of Solution 



For simplicity, we adopt a planar isotropic equation of state, 



(4-3) 



where 7 is the isothermal sound speed squared, as usual. Then, on the disk, equations (2-27), (3-5), and 
(4-2) imply 

Qqv{1 - n)(l + 7) + 7 + 1;^ - n(l + 7) = 0, (4-4) 

"27 + 1 + ^2' 



iV'=-- 



Qo = 2(1 + 7) 



Qo{l + v'') + 2v 



£7 = —2 cot 00. 



(4-5) 

(4-6) 
(4-7) 



We may give equation (4-4) the following quasi-Newtonian interpretation. Because of nonzero pressure, 
the effective gravitational mass-energy is enhanced by the factor (1 + 7) for both the effects of gravitation n 
and the dragging of inertial frames Q[jv{l — n); the net effect of these two terms is balanced per unit inertial 
mass-energy by the "centrifugal term" (with the "radius" scaled out in this self-similar problem) and 
the pressure term 7. Similarly, equation (4-5) is the analog of the Newtonian relationship that the vertical 
gravitational field (cx A^q) just above the surface of the disk is equal to —2ttG times the local surface mass- 
density (oc e); the extra factor represents various relativistic corrections (see the first relation of eq. [4-2]). 
The jump condition (4-6) on the derivative of the frame-dragging term Q'q and the geometrical distortion 
(4-7) of the angular location of the midplane of the disk have, of course, no Newtonian analogs. 

For a cold, slowly rotating, disk, where 7 ^ <C 1, the frame-dragging term oc Q^v is negligible, and 
equation (4-4) recovers the Newtonian approximation for a centrifugally supported Mcstel disk: n w u^. On 
the other hand, from equation (4-10) below, TV' in this limit has a solution which satisfies TV' = at O = 
given by 

2n /I — cos Q ' 



N' « --^ . , (4-8) 

1 n V sm 6 y ^ ' 

with equation (4-7) yielding the location of the disk midplane at 60 ~ 7r/2. Thus, TVq ~ — 2n ~ —2v'^, and 
equation (4-5) now leads to the solution, i w Av^, where E is itself obtained from the first relation of equation 
(3-1) as e « SirTZ^e (see eq. [5-3] below). Thus, we have the identification e « v"^ /2n'R?, which corresponds 
to a Newtonian surface- mass density E = c^TZe/G (radius TZ now having the dimensions of length) related 
to the disk rotational velocity V = cv given by the famous Mestel formula: 



For the fully relativistic situation, the disks are characterized by values of v and 7 that are not very 
small compared to unity. In such an situation, one approach could be to specify these two parameters and 
solve the problem numerically with n and 9o as eigenvalues. In practice, such an approach is very costly. We 
would have to adopt a shooting method in three dimensions for n, Oq and 77 = Q'{0). For nonlinear ODEs, 
the number of operations increases exponentially with the number of eigenvalues. On the other hand, the 
values of v and 7 do not come into play until we get to the disk because the right-hand sides of equations 
(3-10), (3-11), and (3-12) vanish when ^ Qq. Therefore, it is more efficient to treat n and r] as our nominal 
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solution-space parameters, and solve for Qq, v, 7, and e from equations (4-4), (4-5), (4-6), and (4-7) once 
we have integrated the ODEs (3-10), (3-11), and (3-12) for the properties of spacetimc off the disk plane. 

The quantities 80, v 7, and e enter in equations (4-4) -(4-7) with a pattern that allows us to proceed 
as follows. The dynamic equations (3-10), (3-11), (3-12) are cast for 6 ^ 60 as a set of two second-order 
ODEs in N and Q: 



N" = g2 



Q" = Q 



1 — n 
1 + n 



1 — n 
1 + n 



+ {Q'-Q cot e + QN'f - - N' cot 6, 



l + n 



+ {N' - cot Qf -{Q'-Q cot e + QN'f 



- Q'cote. 



(4-10) 



(4-11) 



For given n and ij = Q'(0), we begin with the boundary values, Q = 0, Q' = ?7, iV = 0, A''' = 0, at the pole 
9 = and integrate toward the disk at 6 = ©o (whose value is unknown at this point). At each potential 
choice O for 60, we solve equation (4-4) for 1 -|- 7: 



1-^7 



1 - 7/ 



{Qov + l){l-n)' 

where v is obtained by the following procedure. We first divide equation (4-6) by equation (4-5): 

(1 + -) [Qo(l + r2) + 2r] 



A= — 



N' 



2(l+7)-(l-t;2) 



(4-12) 



(4-13) 



with the value of A known from the off-plane integration to the candidate O for Oq. With the elimination 
of 1 -h 7 through equation (4-12), the last equation implies 



Qov^ + [2 + AQo{l -n)]v + Qo- A{1 + n) = 0. 



This quadratic equation for v yields a solution. 



1 



i+-^(i+„), 

Wo 



(4-14) 



(4-15) 



where we have chosen the sign so that v is well behaved, going to (1 + n)A/2 0, when Qq 0. 
Once V and 7 are known, we can calculate the rescaled energy density 



We have found the disk when @ has a value 60 such that equation (4-7) holds: 



(4-16) 



cot 9o 



(4-17) 



Since the reduced energy density e is positive, we obtain a disk location 7r/2 < 60 < tt. With the other 
obvious limits, 

0<i;<l, 0<7<1, (4-18) 



our parameter space is confined to 



< n < 1, 



r? > 0. 



(4-19) 



5. Results 



Contours of constant j'mn — v space are plotted in Figure 1. Roughly speaking, n is a measure of the 
strength of the gravitational field. It ranges from n = for flat space to n = 1 for maximum rotation. One 
may also argue that since our similarity transformation is r ^ ar, t a'^^^'^H, n cannot exceed unity or the 
passage of time would proceed more quickly deeper in the gravitational potential well, contrary to common 
experience in general relativity. Hence, n lying within the interval (0,1) is an anticipatable result from the 
self-similarity of the basic problem. From Figure 1 we also see that equilibria require lower v for given n if 7 
is larger; this result conforms with the Newtonian intuition that less rotational velocity is needed in a disk 
to offset the self-gravity if there is a greater degree of pressure support. 

For each 7, there is a maximum rotational velocity Vc above which there is no equilibrium. Table 1 
gives the numerical value of Vc as a function of the parameter 7. Our computed value Vc = 0.436 when 
7 = 0.004 is consistent, if we perform a simple extrapolation, with Lynden-Bell and Pineault's (1978b) 
estimate, Vc = 0.438 when 7 = 0, and the latter numbers are what are entered as the first entry of Table 1. 

To understand the last result physically, we examine in Figure 2, for various choices of 7, the behavior 
of the term Qq that governs the dragging of inertial frames in the disk plane. Nonzero values of Qo represent 
how hard one must accelerate to remain at constant (p. In our metric, gu oc (1 — Q"^). Thus, whenever 
Q > 1, an ergorcgion develops and the time-like Killing vector d/dt becomes space- like. Since Q only has 
dependence on 9, the ergoregion is best described by the exterior of a cone whose opening angle 6'orgo is 
defined by Q(0ergo) = 1- Naturally, the ergoregion first appears on the disk when Qo = 1 and ^ergo = ^o- If 
we assume Q is continuous and monotonic, then ^ergo decreases (for the ergoregion above the disk) as Qo 
increases, until finally ^orgo — > as Qo 00. In this limit, when the ergoregion occupies the entire space 
above (and below) the disk, equation (4-4) may be balanced only when n ^ 1, which recovers the upper 
bound on n. The rotation velocity at which Qo diverges is then determined by a limit process on the product 
<3o(l-n). 



5.1. Surface Density of Models 

Table 2 gives the dimensionless coefficient 

«.<i±I*, (5.1, 

(multiplied by 100 to avoid writing too many zeroes) corresponding to a given pair of values v/vc and 7 
that characterizes an equilibrium model. Notice that for small 7 ^ {v/vcY ^ !> we recover the Mestel 
solution, £ w w^/27r. In principle, if one takes the attitude that gas pressures must be three-dimensionally 
isotropic rather than two-dimcnsionally as idealized in this paper, then disks cannot remain vertically thin 
unless 7 <C 1"^. In these restricted circumstances, the entries in Table 2 that violate this constraint are not 
physically self-consistent. Under a broader interpretation of what might be acceptable in the physical world, 
relativistic SIDs might be constructed from non-interacting dark matter particles, in which case allowable 
stress tensors include diagonal forms that are non-isotropic in the sense of equation (2-20). In the Newtonian 
regime, it is also known that strongly magnetized, self-gravitating disks, can be vertically thin even at zero 
rotation speeds (Shu & Li 1997 and references therein). Since there may be useful relativistic analogs of such 
systems, an open attitude would retain all the entries in Table 2 for the sake of mathematical completeness. 

In terms of £, the dimensional surface density of mass-energy £, in units of mass per unit area, as would 
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be measured by a corotating observer, is obtained from 

where TZ is the distance measured from the origin to the point in question along a radial slice at constant cf) 
and t in the equatorial plane ^ = of the disk, 

n = e(^''-^'')/2r. (5-3) 

In equation (5-2), £ is given in Mq km~^ when (? /G is taken to be 0.677 Mq km~^ and TZ is measured in 
km. It should be noted (a) that a proper radius TZ can be defined in the present circumstance because the 

origin neither contains a mass point nor is shielded from an observer at 7?. by a horizon, and (b) that the 
circumference C of a circle in the plane of the disk according to a nonrotating observer at radius TZ is not 
given by 2t:TZ, but by 

C = 2nr = e^^-^'^l'^ 2nn. (5-4) 



5.2. Circular Orbits: Existence and Stability 

We anticipate that many of the more slowly rotating members of the singular isothermal disks (SIDs) 
studied here are unstable to inside-out gravitational collapse in a similar way as their Newtonian counterparts 
(Li & Shu 1997, Shu et al. 2001). We leave for a future endeavor the study of the dynamical, self-gravitating, 

stability of relativistic SIDs and the consequent formation of black holes at their centers if they undergo 
gravitational collapse. Here we ask the simpler question: Are relativistic SIDs kinematically stable in the 
sense of having stable circular orbits for (noninteracting dark-matter) test particles of nonzero rest mass? 
The question is nontrivial, because circular orbits of arbitrary sizes around point masses in the Newtonian 
case are all stable, yet circular orbits lose their stability if they approach too closely the event horizons of 
relativistic point masses (Schwarzschild or Kerr black holes). Is there a similar loss of orbit stability, when 
we go from Newtonian SIDs to relativistic SIDs? 

For a test particle of mass m in the equatorial plane of the disk, the symmetries of the geometry gives 
the conserved quantities: 

Ut = —E = —E/m, u^ = l=l/m. (5-5) 
Since = for an orbit confined to the disk plane, the geodesic equation takes the simple form 

(J) =5'''^[-l-^V* + 2^;V'^-P/1^-2F(r), (5-6) 
where r is the proper time of the particle and V{r) is the effective potential of the problem: 

V{r) = ie-^ne"^" " r-^{Er^-" - Qo^e^"-''")' + Pr'^e^^^-^^"], (5-7) 

and A^O: Po, Qo, and Zq have simple fixed numerical values when N, P, Q, and Z are evaluated in the disk 

plane 9 ^ 9o- 

To have a circular orbit, we need dr/dr = d'^r/dr^ = 0, which implies V{r) — V'{r) — 0. These two 
conditions define the values of specific energy and angular momentum, E and I, needed to yield a circular 
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orbit at radius r. A little algebra shows that the required values of I and E for a circular orbit with r = ro 
are 



where we have defined 



s^o/Vj, (5-8) 



p = 2/(1 -n)-QlT QoVQl + 4n/(l - n)^ 

2(1 -g^) • ^'-'^ 

In equations (5-8) and (5-9), the upper sign choice corresponds to progrado orbits; the lower, to ret- 
rograde ones. For the upper sign choice, the numerator in equation (5-9) goes through zero when the 
denominator does; i.e., for any allowable n, F stays positive when the disk becomes an ergoregion as Qo 
crosses unity. ^Prom our earlier discussions, we recognize the factor as the correct scaling to account 
for the contribution to E of the gravitational "potential energy" per unit mass. We shall prove below that 
F > 1 and that the term \/F represents the contribution of the rest and kinetic energies. Notice now that 
frame dragging adds a positive contribution to the specific energy E of prograde circular orbits, whereas it 
adds a negative contribution for retrograde circular orbits. Moreover, with the lower sign choice in equation 
(5-9), F diverges to +oo when Qo 1—; i.e., the two terms in E cancel to lowest order for large F when the 
disk plane first becomes an ergoregion. Massive test particles in retrograde circular motion have zero total 
energy in this limit. 

The physical interpretation of these results follows from examining the three-velocity of circular orbits 
in the tetrad frame: 

If -I 

where, again, the upper choice corresponds to prograde motion; the lower choice, to retrograde motion. 
Equation (5-10) shows that F is the square of the Lorcntz factor of the particle motion in the tetrad frame: 
F — 1/(1 — Wp). The first term \/F in the expression for E thus represents the usual special-relativistic 
contribution to the rest and kinetic energies of a particle. Moreover, the expression for specific angular 
momentum I in equation (5-8) is now recognized as the velocity Vp in the direction, times the usual 
Lorentz factor correction, times, not the radius TZq, but the circumference of the orbit Cq divided by 27r. The 
above interpretation for F explains why the sequence of retrograde circular orbits terminates when Qo — * 1: 
When the disk develops an ergoregion, massive particles in retrograde motion must travel at the speed of 
light {F — > oo) if they are to resist the dragging of inertial frames. 

No such difficulty affects particles in prograde circular orbit. Equations (5-8) and (5-10) require F > 1 
to make physical sense. We have checked numerically that F as given by equation (5-9) exceeds unity for 
all disks in which the fluid velocity v of the disk is moderately smaller than Vc for any value of 7. When v 
approaches the critical velocity Vc (where the frame-dragging parameter Qo becomes very large), numerical 
errors prevent us from confirming that prograde circular orbits exist. An analytic argument relieves our 
worries on this score. 

For a fixed value of 7, the solution sequence terminates when Qo — > 00. Let us evaluate F in this limit. 
The equation of motion of the disk matter (4-4) may be written as 

l-t>2 

l-n=— — -. (5-11) 

(Qow + l)(l+7) 

Thus, when Qo ^ 00, equation (5-9) with the upper sign choice becomes 

.4 



1 + Vl + 4^2(1 + 7)2/(1 -t;2)2 ^ (5.12) 



which exphcitly satisfies F > 1 for any values of w > and 7 > 0. In the same limit and with the same 
upper sign choice, the three-velocity of the test particle in equation (5-10) is given by 




1-h VI + 41)^(1 +7)V(1-^^^)^ ' 



1+^/1+ 4t;2(l + -y)2/(l 



(5-13) 



Notice that Vp = v when 7 = 0. In other words, the velocity of a test particle in a prograde circular orbit 
equals the velocity of the disk matter when the latter has zero pressure - a satisfying consistency check of 
the result. 

We now wish to investigate whether circular orbits are stable. Imagine perturbing the radial position r 
of the test particle about its equilibrium position vq by a small amount ri , keeping E and I fixed. To lowest 
non- vanishing order upon expansion about ro, equation (5-6) becomes 



Stability of the motion depends on the sign of V"{ro)- If V"{ro) > 0, then circular orbits are stable, because 
there are no perturbations - maintaining the same specific energy and angular momentum - that will produce 
a real solution for ri in the above equation. In this case, of all orbits of a given specific angular momentum 
I, the circular orbit has least specific energy, therefore it is not possible to perturb the circular orbit from its 
equilibrium without giving the particle some additional specific energy, which will then cause it to oscillate 
about the equilibrium radius ro with an "epicyclic" frequency ^/V"{rf)). The radial motion may be pictured 
as rolling up and down the walls of a "valley." On the other hand, if V"{ro) < 0, then circular orbits are 
unstable, because a small perturbation of such an orbit - even one that retains the original specific energy 
and angular momentum - will lead to exponentially growing departures from the equilibrium radius vq. In 
this case, of all orbits of a given specific angular momentum I, the circular orbit has a (local) maximum of 
specific energy, and the test particle becomes unstable by rolling off a "hill." 

In detail, after some algebra, we obtain 



where we have made use of the expressions for I and E from equation (5-8). When Qo = (and hence v — Q), 
1 — n = 1/(1-1-7), P = 1 + 7) and V"{ro) > 0. In this limit, retrograde and prograde circular orbits are 
both stable. In the limit Qq 1 — . F 1/(1 — n^) for prograde orbits, and F{1 — Qq) ^ (1 + r?,)/(l — n) 
for retrograde orbits. Thus, V"{ro) remains positive for both cases. Retrograde and prograde circular orbits 
are still stable forms of motion at the onset of the disk's development of an ergoregion, although even the 
most rapidly moving retrograde particle finds it difficult to resist frame-dragging when Qq ^ As Qo 
passes through unity and approaches 00, retrograde circular motion at velocities less than the speed of light 
becomes impossible, but the product (1 — n)Qo — v"^) / + l)v remains finite and positive, so V"{ro) > 
for prograde circular orbits. We have verified numerically that V"{ro) stays positive between these various 
limits. In summary, prograde circular orbits exist and are stable for power-law isothermal disks from the 
non-relativistic to the ultra-relativistic regime; whereas retrograde circular orbits are possible and stable 
only for disks that do not develop an ergoregion. 

We may state the result in an alternative way that relates to known results concerning circular orbits 
around black holes. Of all orbits of a given specific energy, the circular orbit in a stable/unstable situation 
has a (local) maximum/minimum of specific angular momentum. In the case of a Schwarzschild black hole. 




(5-14) 



V"{ro) = (1 - n)e^'' 



2 n(yF±QoVF^^ +{l-Ql){F-l) 



(5-15) 
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we know that circular orbits that start off at great distances from the event horizon, ro 3> ?'sch) are close to 
the Newtonian limit, and therefore are stable. They remain stable as long as the specific angular momentum 
of the circular orbit continues to decrease with decreasing orbital radius (or circumference). There comes 
a point, ro = 3rsch, when the square of the "epicyclic frequency" F"(ro) changes sign, and the specific 
angular momentum I of the circular orbit has an inflection point, and starts to increase for decreasing tq. 
This violation of "Rayleigh's criterion" signals a transition from stable to unstable circular orbits. Because 
of frame dragging, the case of Kerr black holes is more complicated, but can be similarly elucidated as we 
have done above for the disk case. 

In our power-law disks, every radius ro is similar to any other radius, and spacetime is not flat even 

at infinity. Thus, if a circular orbit is stable at any radius in a given model, circular orbits at all radii arc 
stable. It is hard to imagine how one could realistically construct a self-gravitating disk of rotating material 
otherwise. 



6. Photon Orbits 

The case of photon orbits in the spacetime of relativistic SIDs is also interesting. Our self-similar metric 
equation (2-3) admits a homothetic Killing ^ satisfying 

^i9fi,^ = 2^(/^;i.) = '29^l,^ (6-1) 

In component form, it reads 

+ 9^.C,u + = 25^. (6-2) 

The solution to this equation is 

e = [a-n)t,0,r,0] (6-3) 

Associated with this vector is a conserved quantity F — S.ixk'^ for null geodesies fc** = dx^/dX, where A is an 
affine parameter. Indeed, 

^ = r.^k" = (e^fc'');,fc'^ = = gf^^^k'^k" = ~m\ (6-4) 

where we have used the geodesic equation k^.^k'^ = 0. Therefore, for a massless particle, m = 0, and F is a 
constant of motion. 

In addition, we have the two ordinary Killing vectors associated with the stationarity and axial symmetry 
of the spacetime: 

In total, we have the following three conserved quantities 

E = -kt, l = k4„ F = (1 - n)tkt + rkr ^ kr = ^[T + {1 - n)tE] (6-6) 

The null condition k^k^ = can be used to determine ke: 

klg'' + 2ktk^g"' + k^'^ + klg"-^ + klg"' = 0, (6-7) 

which implies 

ke = ±-i= + 2Elg'^ - l^g^l' - 1 [F + (1 - n)tEf ff-j 



1/2 
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= ±e^/2 { [Er^-^e-'' - g/e"^]' - /^e'^^ - [F + (1 - n)Etf e"^}'^' . 
Finally, the geodesic is described by 

= r-^^+^h-^'EQ + r-2Ze^-2^(l - Q^), 
j^r ^ ^N-z^-i [r + (1 - n)Et] , 

= ±r-2e^-^/2 { [£ri-"e-^ - QZe"^] ' - Z^e'^^ - [F + (1 - n)i;t]' e"^} 
Divide everything by fc*, and we have 

d4> _ Q + a{l - Q^) ^^_,^_p 



1/2 



dt 



1-Qa 



-r e 



d\nr e^N-z^n-ip 



dt l-Qa ' 

6^-^/2(1 + n)r"-i 



[(l-gaf-a2-/32e 



2N-Z 



1/2 



where 



+ (1 - n)t 



„n-l 



To make this system autonomous, we extend the space to include a and /3 as two of the variables 
self-similarity of the problem makes it convenient to define the reduced radius and time as 



We then have 



C = lnr and dr = r"-~^ dt. 

# _ Q + a{l- Q^) j^_p 
dr 1 — Qa ' 

dC _ 

dr 1 — Qa ' 

N-Z/2 



— = ± 1 + n — 

dr 1 — 



rf7 



= (l-n) 



[(1 - - - P'e 

2N-Z 02 



2 „2iV-Zl 1/2 



1- 



^2N-Z 



1-Qa 



1-Qa 

N-Z/2 



2 2N-Z-1 1/2 



1-Qa 



[(1 _ Qaf -a'- /?2e2^-^] 



dt 
d^ 



=,f(i-») 



(6-8) 

(6-9) 

6-10) 
6-11) 

6-12) 

6-13) 

6-14) 
6-15) 

6-16) 

, The 

6-17) 

6-18) 
6-19) 
6-20) 
6-21) 
6-22) 
6-23) 
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6.1. Dynamics and Geometry of the Photon Trajectories 

We need to specify initial values for the six dependent variables (</>, (, O, /3, a, t) and integrate forward in 
T. Without loss of generality, stationarity, axial symmetry, and self-similarity imply that we may take t = 0, 
= 0, and r = \ (or ^ = 0) at r = 0. We shall also assume that all photons begin by being emitted from the 
plane of the disk, 8 = 8o- Except for a set of measure zero (involving photons travelling outward exactly 
along the rotation axis), we shall find that this last assumption also results in no loss of generality because 
even photons emitted by external sources outside of the disk are soon bent to cross the disk plane. This 
behavior can be attributed to the existence of an adiabatic invariant J (see below), which places a constraint 
on the trajectories different from those presented by the classical integrals I, E, and V. In any case, we are 
now left with only two arbitrary specifications, the initial values a = a* and /? = /3,, from which equations 
(6-16) allow us to reconstruct the two constants of motion, 

l/E = a^, T/E = f3^. (6-24) 



Only the combinations l/E and T/E provide isolating integrals because the principle of equivalence 
forbids photons having values of /, E, and T that differ by only a single multiplicative factor from having 

fundamentally different trajectories in spacetime. By this method of counting, wc sec that (self-similar) 
photon trajectories are completely determined by well-behaved integrals of motion, and thus chaos does not 
enter the problem. Indeed, self-similarity and axial symmetry decouple the three ODES for G, a, and f3 from 
the rest of the set of equations (6-18) (6-23), so that the other three variables {4>,(,t) may be computed 
by post-processing. Although useful for proving the absence of chaotic photon orbits in this problem, the 
above discussion unfortunately lends little descriptive power to the geometry of the photon orbits. Instead, 
we adopt the following approach for choosing the parameters of the initial conditions. 

Let us consider the photon orbit as seen by a locally nonrotating observer (LNRO) whose basis one-forms 
are given by 



,(*) 



a;(*)=re^-~/2(d,^-r" 



Qdt), 



(6-25) 
(6-26) 
(6-27) 
(6-28) 



Projected onto this frame, fc^^' — k • w^^'. In component form, 

= [Er'-^e-'' - e^^^lQ] = Er-^er^'/^l 

fc(r) ^ e(JV-2)/2^-i[r + (1 _ n)Et] = i?e(^-^)/2r-"/3, 
kie) = ±r-ieiv/2 I [Er^-n^-N _ g^^-P] 2 _ ^2g-2P _ [r + (1 . 

r 2 , /o2 2JV-Z-| 1/2 

Notice that since A:*^*' is the energy measure by the LNRO, it is always positive. That means E and 1 — Qa 
have the same sign and we can safely omit the absolute sign in k'^^K Let {tp, x) be the direction of a photon 
trajectory seen by this LNRO. Then these angles are related to equations (6-18)-(6-23) by 



Qa), 



n)Etf e^-^ 



1/2 



(6-29) 

(6-30) 
(6-31) 



(6-32) 



cosV; = fc(^V^ 



(t) 



± 



a 



2^2N-Z 



(1 - Qaf 



1/2 



(6-33) 
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tanx = fc(^Vfc^'' = (6-34) 

It is more intuitive to specify tp and x as initial conditions. In particular, if we adopt the other choices 
discussed at the beginning of this subsection, we have at r = 0: 

= 0, C = 0, e = eo, (6-35) 

P = e^^^--^^^^^p^^P., a=-^^^p^^a., t = 0. (6-36) 
y sin i/» sin X ± 1 Q sin i/» sm x ± 1 



6.2. Sign Choices 



To solve the geodesic equations (6-18)-(6-23) numerically, we need to be careful about the ± signs and 
the term in the square root of equation (6-20), which we call A: 



A s (1 - Qay -a^ - j3^e 



2 2N-Z 



(6-37) 



In order for dQ/dr to be real, we need dA/dr to vanish whenever A attains zero from positive values 
(otherwise A can become negative). By straightforward differentiation, it is easy to show 



2 dr dr dr 



(6-38) 



where we have used that dO/dr = when A = to eliminate derivatives of functions of only 8. Upon 
substitution of equations (6-21) and (6-22), the last expression becomes 



[Q + a{l - Q')] a{n - 1)^^ + Pe'''-^{1 - n) 
'- ■' 1 — Ua 



O'c 



1 - 



1-Qa 



With A = in equation (6-37), the term in the last square bracket equals 

"(1 - Qa)Qa + q^' 
1-Qa 



(6-39) 



(6-40) 



and when this is substituted into equation (6-39), we find that the first and second terms algebraically cancel. 
Thus, dA/dr vanishes when A = 0, which is the desired result. 



Initially, we require that for ^ < tt/2, dQ/dr is negative: 

de 



dr 



= -{l + n)e^-^/^cosi) 



Similarly, we want dQ/dr positive whenever — 7r/2 < x < t^/"^- This requires, after a little algebra, 

^ = ie^-^/^sinV'cosx > 0, 
dr 



(6-41) 



(6-42) 



i.e. we choose the -|- sign in the initial conditions for a and /3. We can then integrate 6 forward in r until 
A = 0, where we change the sign of dQ/dr (and correspondingly the sign in da/dr). 



- 18 - 



6.3. Results 

First we shall discuss the photon orbits in a spacetime without the ergoregion (Qo < !)• As an example, 
we choose the parameters n = 0.4 and 7 = 0.5 (or equivalently v = 0.16), and consider photons launched in 
different directions as seen by the LNRO. We assume that photons cross the disk plane without absorption 
or scattering, an assumption that is more likely to apply to massless low-energy neutrinos than real photons, 
for which a disk with surface density 10^^ g cm~^ (if we are talking about stellar-mass disks of sizes ~ 1 
km) is not likely to be optically thin. Our formal usage of the phrase "photon orbits" must henceforth be 
understood either as a metamorphic rather than literal device, or as applicable only while the photon is 
travelling above or below the disk plane. 

We find that outgoing photons with — 7r/2 < x < 1^ spiral out to ^ = 00, while ingoing photons with 
X > 7r/2 or X < — 7r/2 reach a minimum radius and then spiral out to C = 00 again. Only one photon 
launched along the disk with i\) = 7r/2) and a = reaches the origin. This was first discovered by Lynden- 
Bell & Pineault (1978b) for the cold disk. Photons launched in the retrograde direction, reach a minimum 
coordinate angle (jimm and are then dragged to go in the forward direction. Figure 3 shows these orbits. 

We now discuss the out-of-plane behavior of the orbits. No matter what is the initial condition, almost 
all orbits are focused and eventually trapped by the disk. A (noninteracting) photon will typically penetrate 
the disk many times. Each time it reaches a turning point &i, where i labels the number of penetration (see 
Fig. 6). The only trajectory that can escape falling into the disk is one launched exactly along the rotation 
axis. However, such a photon suffers an infinite redshift as it propagates away from the disk, and an observer 
located at an infinite distance above or below the disk cannot see it. (This result was first discovered for 
counter-rotating disks, and is probably generic to self-similar configurations that contain an infinite total 
mass.) 

When interactions with the disk are ignored, it appears that Oj decreases with each disk crossing i. The 
result can be demonstrated to arise from adiabatic invariance. The relevant invariant is easily computed 
conceptually. The conjugate momentum to 9 is kg. Hence we can define the action integral, 

J = jkede^j^\E\ y°""%i-"{[l-0a]'e^-2iV-a2e2-2Ar_^2|i/2^Q_ (g_43) 

Here we may treat a and /? as functions of O along the photon's trajectory during the current cycle. The 
quantity r{Q) varies slowly over this one cycle because it is a monotonic function which does not oscillate. 
Hence, we may approximate it by its mean value. To zeroth order, we may take r(6) outside of the integral. 
Thus 

J = |£;|ri-" [J(e„,ax) - J(e„,in)] . (6-44) 

where 

J(e„0 ^ |[1 - Qaf e^-^" - a^e^-^^ - fi'V" dS. (6-45) 

1 + n 7o L J 

Because the integrand is a known (numerical) function of O, we can tabulate the integral as a function of 

its argument for each photon trajectory. We can then invert the expression equation (6-44) to recover 

C = Inr as a function of 6i„i„, and ©max- 

C = In [J(©max) - J(©min)] + C, (6-46) 

1 — n 

where C is the constant \n{J/E) divided by 1 — n. The resulting curve is plotted as a dashed locus in Figure 
6. The concordance between the dashed curve and the actual envelope of the solid photon trajectory is a 
measure of the goodness of the adiabatic invariant J. 
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It might be argued that the vahdity of the focusing effect on photons seen in Figure 6 is compromised by 
the assumption that photons continue in their original direction when they cross the disk plane. However, 
since all orbits, except for a set of measure zero, always return to the disk plane, the qualitative effect 
remains the same even if we were to allow photons to interact with the disk matter. The absorption and 
rc-cmission or the scattering of photons when they cross the disk would result in a slow transfer of photons 
from the inner disk to the outer disk because no photon can permanently escape from the disk. The focusing 
of photons toward the equatorial plane is perhaps a generic feature of relativistic disks (and not just peculiar 
to these self-similar configurations), and may present an obstacle to some classes of beamed-jet models for 
gamma-ray burst sources. 

Photon orbits in a spacetime with Qo > 1 are more complicated than those discussed so far. As an 
example, let us take n = 0.75, v = 0.32 (and 7 — 0.5). All forwardly propagating photons (x > 0) escape 
to C = 00. The ones launched outward (x < tt/2) escape directly, with the 6 dependence mimicking the 
behavior of those in spacetime without an ergoregion. The ones launched inward (x > 7r/2) reach a minimum 
radius and then escape. In the interior region, the photons spiral toward the axis above and below the disk. 
The turning points in 6 increase as the photons approach the minimum radius. Once past the starting point 
= Inr = 0, we see the familiar 6 behavior governed by the adiabatic invariance described above. Figure 7 
shows the C and 6 behavior of these trajectories. 

The backwards photons are divided into two classes, separately by the surface Xc(V')- Roughly speaking, 

this surface corresponds to _E = in the outgoing direction. All ingoing photons fall towards the origin 
directly. For an outgoing photon, if x > Xc) then it escapes to infinity. On the other hand if x < Xc) the 
photon reaches a maximum radius, and then falls to the origin (see Fig. 8). For the spacetime consider 
here, the E = surface is plotted in Figure 9. Lynden-Bell & Pineault (1978b) gave an analytic expression 
for Xc(7r/2) = sin^^ (— l/Qo)- For this particular spacetime, Qo ~ 4.41. Therefore Xc(7r/2) = —0.073, which 
agrees with our numerical result. 

7. Conclusion 

We have solved by semi-analytic means the Einstein field equations for axisymmetric, self-similar, rel- 
ativistic disks with "flat" rotation curves, including finite levels of pressure support. These spacetimcs are 
not asymptotically flat and cannot describe correctly the behavior of isolated astrophysical objects when 
examined at distances that are very large compared to their natural gravitational radii. Nevertheless, the 
solutions may yield some insight into the near-field solutions of rapidly rotating, compact objects. 

As expected from first principles, the solution space is parameterized by two dimensionless numbers, 

V and 7, that represent the disk rotation speed and the square of the isothermal soimd speed when both 
are normalized appropriately by the speed of light c. The qualitative behavior of these disks resemble those 
found by Lynden-Bell & Pineault (1978b) for the cold disk. This is encouraging because cold disks are known 
in their Newtonian limits to be violently unstable to a wide variety of spiral and barlike perturbations, and 
we cannot expect their relativistic counterparts to behave much better. A proper stability analysis of the 
disks discussed in the current paper remains a task for the future. 

Ergoregions develop for relatively low rotation velocities in our disks and take the shape of (the outside 
of) a cone centered around the axis. As the rotational velocity increases, the "ergo-cone" closes up towards 
the axis. For each 7, there is a maximum velocity Vc beyond which no equilibrium can exist due to infinite 
frame dragging. It should be noted that this maximum velocity lies well below the special relativistic limit 



- 20 - 



of the speed of light. 

We examined the behavior of test particles with nonzero rest mass placed in circiilar orbits in the plane 
of the disk. We found that prograde circular orbits exist and are stable for the full range of disk models 
in this paper. Retrograde circular orbits are also stable when they exist, but cannot be maintained against 
frame dragging by particle velocities less than the speed of light when the disks develop ergoregions. 

We also carried out a systematic study of planar and non-planar photon orbits. Most interestingly, we 
found that all photon orbits are ultimately attracted toward the plane of the disk because of the operation 
of a general adiabatic invariant. Although the formal result depends on the disk being optically thin to 
the propagating photons (an unlikely state of affairs), we gave physical arguments why the generic effect 
may pose defocusing difficulties for some classes of models of gamma-ray burst sources that rely on beamed 
jets along the rotation axis of rapidly rotating compact objects. To be sure, the effect in realistic flattened 
systems that do not have infinite mass and spatial extent may be less dramatic than the one found here 
for relativistic SIDs. A lower bound on the effect might be obtained by examining the analogous properties 
of photon orbits in a Kerr geometry, which in other respects mimics the spacetime analyzed in the current 
paper. 

It is our belief that the current investigation has just begun to scratch the surface of a potentially 
very rich mine for general relativity to explore. As discussed in the Introduction, the study of self-similar 
(Mestel) disks in the Newtonian limit has uncovered rich veins relating to the stability and collapse of such 
objects that have illuminated astronomers' understanding of real-world objects such as protoplanetary disks 
and triaxial and spiral galaxies. In addition to serving as useful testbcds for numerical relativity codes, the 
relativistic generalization of such studies could shed light on topics such as the efficiency of gravitational 
radiation and the possible generation of naked singularities during gravitational collapse. 

We thank Donald Lynden-Bell for useful conversations. This work has been supported by the National 
Science Foundation by a Graduate Fellowship awarded to MJC and by grant AST-9618491 awarded to FHS. 
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Fig. 1. — The gravity index n vs rotational velocity v for different values of the pressure parameter 7. The 
dashed line is the empirical approximation to n for 7 = 0, found by Lynden-Bell and Pineault. 
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Fig. 2. — The frame dragging parameter Qo vs rotational velocity v for different values of the pressure 
parameter 7. 



Table 1. Critical Velocity as a Function of Sound Speed Squared 
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(a) Azimuthal angle '/>(t) for photons launched in different directions 




-40 I , , , \ , , , \ , , , I 

20 40 60 

T 

(b) Logarithmic radius C(r) for ingoing and outgoing photons 



Fig. 3. — Behavior of photon trajectories without the ergoregion. 



Fig. 4. — A typical photon trajectory. Here the distance is plotted on a log scale, normalized such that the 
closest approach to the origin is 1. The nonoscillatory curve at the bottom is the projection of the photon 
trajectory onto the equatorial plane of the disk. 



Tabic 2. 100 f as a Function of v jvc and 7 
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Fig. 5. — A photon launched along G = 10 ^ is being trapped by the disk. Again, the distance is plotted 
on a normalized log scale. 
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Fig. 6. — Photon orbit showing as a function of C. The turning point in Q decreases each time after the 
photon penetrates the disk. The sohd Unc is from solving geodesic equation directly, while the dotted line is 
the computed envelope from adiabatic invariance. 
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(a) Logarithmic radius (^{t) for ingoing and outgoing photons 




10 20 30 40 50 60 

(b) Normahzed polar angle 0*(t) for ingoing photons. Here 0, is normalized such that the disk is located 

at 7r/2. 

Fig. 7. — Behavior of photon trajectories lauched in the forward direction with an ergoregion. 
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Fig. 8. — Logarithmic radius C(t) for photons on either side of i? = surface 
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